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We present results for the nucleon axial charge qa at a fixed lattice spacing of l/o = 1.73(3) GeV 
using 2+1 flavors of domain wall fermions on size 16 3 x 32 and 24 3 x 64 lattices (L = 1.8 and 2.7 fm) 
with length 16 in the fifth dimension. The length of the Monte Carlo trajectory at the lightest m n is 
7360 units, including 900 for thermalization. We find finite volume effects are larger than the pion 
mass dependence at = 330 MeV. We also find that gA exhibits a scaling with the single variable 
m w L which can also be seen in previous two-flavor domain wall and Wilson fermion calculations. 
Using this scaling to eliminate the finite-volume effect, we obtain gA = 1.20(6)(4) at the physical 
pion mass, m-K = 135 MeV, where the first and second errors are statistical and systematic. The 
observed finite-volume scaling also appears in similar quenched simulations, but disappear when 
V > (2.4 fm) 3 . We argue this is a dynamical quark effect. 

PACS numbers: 11.15.Ha, 11.30.Rd, 12. 38. Aw, 12.38.-t 12.38.Gc 



The isovector axial charge of nucleon is a fundamen- 
tal observable in hadron physics. It is defined as the 
axial vector form factor at zero four-momentum trans- 
fer, g A — Ga(0). The axial vector form factor is 
given by the nucleon matrix element of the axial vector 
current, A" = vjj ^75(1"° / , with u,d quark doublet 
if), (n>\A«\n) =u n ,l^G A (q 2 ) +iq^G P (q 2 )h 5 (T a /2)u n , 
where Gp is the induced pseudoscalar form factor, r a 
an isospin Pauli matrix, and q the momentum transfer, 
q^ = P^—p 1 ^ ■ Experimentaly g A has been obtained very 
precisely, gA — 1.2695(29) through neutron @ decay. 

g A is related to the spontaneous breaking of the chi- 
ral symmetry of the strong interaction through the well- 
known Goldberger-Treiman relation This relation 
shows that gA is proportional to the strong pion-nucleon 
coupling at q 2 « 0. Furthermore the Adler-Weisberger 
sum rule 0, |3| reveals that g A differs from unity for a 
structure-less nucleon through the difference between the 
integrals of the total cross sections of the n + p and n~p 
channels. These are in good agreement with experiments. 

Hence the axial charge allows us to perform a preci- 
sion test of (lattice) QCD in the baryon sector. More- 
over, since it is an isovector matrix element, only con- 
nected quark diagrams at q 2 = contribute, making 
the calculation technically simpler. Furthermore, the 
renormalization of the axial vector current is simplified 
when utilizing a lattice chiral fermion action such as 
the domain wall fermion (DWF) action 0, d, Q- How- 
ever, it is difficult to control finite- volume effects (FVE's) 
in gA, as suggested by quenched calculations [8]. The 



FVE's have been investigated in effective models [9T [Io|. 
and also in heavy baryon chiral perturbation theory 

(HBChPT) [HEEm. The FVE ' S in HBChPT are in- 
consistent with lattice calculations unless contributions 
of the A baryon resonance are included. Such con- 
tributions introduce several unknown parameters, and 
the FVE seems sensitive to one of them, the A-iV cou- 
pling ll|, [l3j, [l4j. Until now there has been only one 



investigation in the light quark mass region with 2+1 fla- 
vors, the mixed action calculation [l5| by LHPC which 
reported the FVE is small for a volume of V = (2.5 fm) 3 . 

In this letter we report a result for g A using 2+1 fla- 
vors of dynamical DWF with several quark masses corre- 
sponding to m v = 0.33-0.67 GeV, on two volumes with 
spatial size L — 1.8 and 2.7 fm. Our light quark masses 
and the spatial volumes allow a detailed investigation of 
the FVE in full QCD where the sea and valence quarks 
are identical. It is found that the FVE is not negligible at 
the lightest m„, even on the larger volume, and that g A 
exhibits a scaling behavior with the single variable m n L. 

Our calculation is performed with a fixed lattice spac- 
ing, two space-time lattice sizes, 16 3 x 32 and 24 3 x 64, the 
Iwasaki gauge action [l6| with f3 — 2.13, and the DWF 
action with a fifth dimension of size 16 and M5 = 1.8. 
Each ensemble of configurations uses the same dynam- 
ical strange quark mass, m^ ca = 0.04 in lattice units. 
We use four light sea quark masses on the 24 3 lattice, 
TOsca = 0.005,0.01,0.02, and 0.03 and the same heav- 
iest three masses on the 16 3 lattice. The ensembles, 
described in 17, UH, were generated using the RHMC 
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TO/ 



0.005 



0.01 



0.02 



0.03 



TO w [GeV][21] 0.3313(13) 0.4189(13) 0.5572(5) 0.6721(6) 
(2.7 fm) 3 1.083(50) 1.186(36) 1.173(36) 1.197(30) 
(1.8 fm) 3 N/A 1.066(72) 1.115(58) 1.149(32) 

TABLE I: g A and (V = (2.7 fm) 3 only). 

algorithm [19j | with trajectories of unit length. The mea- 
surements were performed at the unitary points only, 
mf = m V ai = nisca- We use the mass of the £l~ baryon to 
determine the inverse of the lattice spacing 1/a = 1.73(3) 
GeV 2^, 21 1 . The residual quark mass due to the finite 
size of the fifth dimension is 0.00315(2). The non-zero 
lattice spacing error is small in our calculation because 
the DWF action is automatically off-shell 0(a) improved. 

Four measurements are carried out for the 24 3 ensem- 
bles on each configuration. The number of Monte Carlo 
trajectories used for measurements is 6460, 3560, 2000, 
and 2120 for m f = 0.005, 0.01, 0.02, 0.03, respectively 
with 10 trajectory separations for to/ — 0.005, 0.01 and 
20 for 0.02, 0.03. The measurements are blocked into 
bins of 40 trajectories each to reduce auto-correlations. 
On the 16 3 ensembles we use 3500 trajectories separated 
by 10 trajectories at to/ = 0.01, and 0.03, and by 5 at 
0.02. The data are blocked with 20 trajectories per bin. 

The axial charge is calculated from the ratio of the ma- 
trix elements of the spatial component of the axial vector 
current and the temporal component of the vector cur- 
rent, V t a = ^7t(r a /2)V', (n'\A^\n)/(n'\V t a \n) = g A - This 
ratio gives the renormalized axial charge because and 
Vp, share a common renormalization constant due to the 
chiral symmetry of DWF. In our simulation the two con- 
stants are consistent to less than 0.5% at the chiral limit. 
In order to increase the overlap with the ground state, 
the quark propagators are calculated with gauge invari- 
ant Gaussian smearing [22j and we employ sufficient sep- 
aration in Euclidean time, more than 1.37 fm, which 
is the lar gest used so far in dynamical calculations of 
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23l . |24| , between the location of the nucleon source 
and sink to minimize excited state contamination. 

The plateaus of gA computed on volume V = (2.7 fm) 3 
are shown in Fig. [T] We checked that consistent results 
are obtained by either fitting or averaging over appro- 
priate time slices, t = 4-8, and also by fitting the data 
symmetrized about t = 6. The larger volume data can 
be symmetrized because the source and sink operators 
are identical in the limit of large statistics. We note that 
the length of our lightest mass run is already the longest 
we know of for comparable simulation parameters. Re- 
sults obtained from the fit using the unsymmetrized data, 
presented in the figure with one standard deviation, are 
employed in the analysis. 

Figure [2] shows our result for gA- The results are also 
presented in TableQ] The (2.7 fm) 3 data are almost inde- 
pendent of the pion mass (squared) except for the lightest 
point which is about 9% smaller than the others. A set 
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FIG. 1: Plateaus of g A . V = (2.7 fm) 3 and to/ = 0.005, 0.01, 
0.02, and 0.03, from top to bottom. 
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FIG. 2: gA- Dashed and solid lines denote the fit results and 
chiral extrapolation in infinite volume, respectively. The open 
circle is extrapolated result at m n = 135 MeV. 



of the results obtained with a smaller volume, (1.8 fm) 3 
shows a similar downward behavior, albeit with relatively 
larger statistical uncertainties. An earlier two flavor cal- 
culation by RBC FJ with spatial volume (1.9 fm) 3 and 
1/a = 1.7 GeV showed a clear downward behavior, but 
it sets in at heavier pion mass. 

We suspect that this pion mass dependence driving 
gA away from the experimental value is caused by the 
finite volume of our calculation: in general such an ef- 
fect is expected to grow as the quark mass gets lighter 
at fixed volume, or the volume decreases for fixed quark 
mass. More quantitatively we observe in the figure that 
the two flavor result with V — (1.9 fm) 3 significantly de- 
creases at m 2 « 0.24 GeV 2 , while the 2+1 flavor results 
with V = (2.7 fm) 3 do not decrease even at to 2 « 0.17 
GeV 2 . Another trend of the FVE seen in Fig. [5] is that 
all the 2+1 flavor, smaller volume data are systematically 
lower than the larger volume data. Similar behavior was 
observed in quenched DWF studies [8j, l25[. However, for 
pion masses close to our lightest point such a sizable shift 
is not observed when V is larger than about (2.4 fm) 3 , 
not only in the quenched case, but also the 2+1 flavor, 
mixed action, calculation in [15| . as shown in Fig.[2j On 
the other hand, our results suggest that V = (2.7 fm) 3 
is not enough to avoid a significant FVE on gA when 
TOtt < 0.33 GeV in dynamical fermion calculations. 

In order to more directly compare the various results, 
we plot gA against a dimensionless quantity m^L, in the 
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FIG. 3: m^L scaling of gA- Top, middle, and bottom panels 
are dynamical and mixed action DWF, dynamical Wilson, 
and quenched DWF results, respectively. In bottom panel, 
open symbol is same as in top panel. 

top panel of Fig. [3] Interestingly, we find that the 2+1 
flavor results on both volumes and the two flavor ones 
reasonably collapse onto a single curve that monoton- 
ically increases with m^L; in other words, they exhibit 
scaling in m^L. A similar scaling also appears in dynam- 
ical two flavor (improved) Wilson fermion calculations as 



shown in the middle panel of Fig. [3] 13j, [23J, |24[ for the 
unitary points K sca = K va i, with various volumes (0.95- 
2.0 fm) 3 , pion masses 0.38-1.18 GeV, and gauge cou- 
plings. The large difference in Wilson data at j3 = 5.60 
on m % L ~ 6.5 can be described by different choices of the 
renormalization constant of A^. While the trend is simi- 
lar in the quenched DWF case H, [25| with pion masses in 
the range 0.39-0.86 GeV and l/a = 1.3 GeV (see bottom 
panel, Fig. [3]), the scaling is violated for the point with 
smallest m^L on V = (2.4 fm) 3 . The lightest point does 
not follow the (1.8 fm) 3 data: they differ by 2.5 standard 
deviations (er) at m^L ~ 5, suggesting that there are 
non-universal terms that depend separately on and 
V. In particular, this effect may be due to the presence 
of a quenched chiral log [SiJ . From [26] , the size of the 
effect at this mass can readily explain the discrepancy 
observed with the dynamical m n L scaling. Note, at this 
mass, but going to V — (3.6 fm) 3 , no FVE is detected in 
the quenched case as can be seen in Fig. [2j 

It is interesting to compare our larger volume, 2+1 fla- 
vor result with the mixed action, 2+1 flavor result with a 
similar volume [HI , denoted by the left triangle in the top 
panel of Fig. [3) At heavy pion mass the results are sta- 
tistically consistent and essentially independent of m n L. 



At m^L ~ 4.5 the mixed action result, however, is larger 
than ours by (a combined) 2.1<r, but consistent with the 
quenched DWF result with (2.4 fm) 3 volume Q, (the up 
triangle in the figure), which is also 2.2er larger than our 
result. The mixed action result was also calculated with 
an even larger volume, and no FVE was detected [l5j |. 
Again, similar to the quenched case. 

A possible explanation of the difference is that it is sim- 
ply a dynamical fermion effect. The mixed action results 
are partially quenched, calculated with improved stag- 
gered sea quarks and domain wall valence quarks. The 
pion mass of valence DWF is tuned to match the lightest 
pseudoscalar meson mass computed entirely with stag- 
gered fermions. However, there is an ambiguity 27], l28[ 
in choosing the staggered pseudoscalar meson for this 
tuning since there are several. Thus, if the valence quark 
is effectively much lighter than the sea quark, a mixed 
action calculation effectively becomes quenched. This 
may lead to a deviation from the unitary calculation and 
consistency with the quenched calculation. Mixed ac- 
tion ChPT shows the presence of partially-quenched log s 
whose size is consistent with the observed effect 29|, 30j . 

For the chiral extrapolation of gA, we attempt to in- 
clude the FVE in our data. While the pion mass depen- 
dence of gA, including the FVE, has been investigated in 
the small scale expansion (SSE) scheme of HBChPT [l3j ]. 
the size of the FVE on V = (2.7 fm) 3 predicted in SSE is 
less than 1% in our pion mass region. The correction is 
much too small to account for the observed FVE in our 
data. This suggests that the FVE in HBChPT, which is 
estimated by replacing all loop integrals by summations, 
is not the leading FVE in gA, as one in e regime [3l| . 
We also note that our attempts to fit the mass depen- 
dence of the data to HBChPT failed, which is likely due 
to the heavier quark mass points being beyond the radius 
of convergence of ChPT HE Is! • 

Instead of the SSE formula, we assume the following 
simple fit form, including the FVE in a way that respects 
the scaling observed in the data, A + BmJ + C/y (m^-L), 
with fv{x) = e~ x , and where A, B, and C are fit pa- 
rameters. We employ a constant term and one linear 
in the pion mass squared, the leading contributions to 
gA in infinite volume. The third term corresponds to the 
observed FVE, taken as a function of m^L only, and van- 
ishes rapidly towards the infinite volume limit, L — > oo 
at fixed pion mass. The same m^L dependence appears 
in one of the FVE contributions in Ref. . We note that 
this simple form is used to estimate the FVE's in the data 
but not the value of gA in the chiral limit at fixed L. In 
the end, we chose this simplest form, in part, because 
the fit result at the physical point is not sensitive to the 
particular choice of fv(x), as discussed below. 

In Fig. [2] we see that the 2+1 flavor data are described 
very well by this simple fit (x 2 /d.o.f. = 0.45), using 
data computed on both volumes simultaneously. The 
L — > co extrapolation (solid line) in turn allows an ex- 
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trapolation to the physical pion mass (m^ = 135 MeV), 
c]A = 1.20(6)(4), where the first error is statistical. The 
second error is an estimate of the systematic error de- 
termined by comparing this result with that from fits 
using different choices of fv(x), e.g., the full form in 
x~ 3 , and ml e~ x /x 1 / 2 . The latter is similar to HBChPT 
when m^L ^> 1 [ill [l2l Il3j |. The extrapolated value 
is not sensitive to the choice of fv, and is also consis- 
tent with a linear fit to the three heaviest points on the 
larger volume, — 1.17(6). We also fit our data, with 
and without the fv term, to the 2-loop formula from 
HBChPT [32} and find that the extrapolated result is 
less than 1 and that the fits are generally unstable. This 
is because the many unknown low energy constants can- 
not be determined well from only four data points, even if 
some of them arc fixed. More importantly, though the 2- 
loop formula extends the range of the chiral expansion, it 
is still only large enou gh t o include our lightest point, as 
demonstrated in Ref. [32j. The systematic error arising 
from the difference of the renormalization constants for 
and U M is much smaller than the quoted systematic 
error. From the fit result with fv(x) = e~ x , we estimate 
that if one aims to keep FVE's at or below 1%, then for 
TOtt = 330 MeV, spatial sizes of 3.4-4.1 fm (m^L w 5.7- 
6.9) are necessary. 

We have computed the nucleon axial charge with 2+1 
flavors of DWF on a largest volume of (2.7 fm) 3 , a lightest 
TOtt of 330 MeV, and large statistics. We have observed 
a downward dependence of gA on our larger volume, 
and by comparing our results with those using differ- 
ent volumes, numbers of flavors, and lattice fermions 
as a function of the single variable m^L, concluded it 
is caused by the finite volume used in our calculation. 
The data appear to scale well in this variable except in 
the quenched and partially quenched, mixed action cases 
when the volume exceeds roughly (2.4 fm) 3 , which may 
be due to the presence of (partially-) quenched chiral logs. 
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